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ON A DEGENERATE NON-LOCAL PARABOLIC PROBLEM DESCRIBING 
INFINITE DIMENSIONAL REPLICATOR DYNAMICS 


NIKOS I. KAVALLARIS, JOHANNES LANKEIT, AND MICHAEL WINKLER 

Abstract. We establish the existence of locally positive weak solutions to the homogeneous Dirichlet 
problem for 

ut = uAu + u |Vu|^ 

Jn 

in bounded domains D C R"' and prove that solutions converge to 0 if the initial mass is small, whereas 
they undergo blow-up in hnite time if the initial mass is large. We show that in this case the blow-up 
set coincides with D, i.e. the hnite-time blow-up is global. 
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1. Introduction 

In a bounded domain Q. C N > 1, we consider nonnegative solutions to the quasilinear degenerate 
and nonlocal parabolic problem 

ut = uAu + u j |Vup, xEH , t > 0, 

^ u{x,t)=i)^ X G (90, t > 0, (1-1) 

, u{x,o) = uq{x), X g n, 

which arises in a game theoretical description of replicator dynamics in the case of a Bomze-type 
inhnite dimensional setting [8] by pursuing a modeling procedure introduced in [231 [Ml [38] and which 
actually assumes steep payoff-kernels of Gaussian type. For completeness in this direction we include 
a concise derivation of the particular parabolic equation in (II.ip in the Appendix A. 

Strongly degenerate diffusion meets nonlocal gradient sources. From a mathematical 
perspective, the evolution in (|l.ll) is governed by two characteristic mechanisms, each of which already 
gives rise to considerable challenges on its own. Firstly, diffusion in (II.ip is strongly degenerate at 
small densities in the sense that near points where u = 0 typical diffusive effects are substantially 
inhibited. Indeed, already in the unforced counterpart of (jl.ll) with general power-type degeneracy, 
as given by 

Uf = vPAu (1-2) 

with p > 0, it is known that the particular value p = 1, corresponding to the choice in (|l.ll) . marks 
a borderline between somewhat mild degeneracies and strongly degenerate diffusion: In the case 
when p < 1, namely, (|1.2p allows for a transformation into the porous medium equation vt = Ax™ 
with m := > 1, thus meaning that in this case unique global continuous weak solutions to the 
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associated Dirichlet problem exist for all reasonably regular nonnegative initial data ([3]), and that 
these eventually become positive and smooth, and hence classical, inside ([6]). If p > 1 then 
nonnegative global weak solutions can still be constructed for any nonnegative continuous initial data, 
but they need no longer be continuous ([5]) nor uniquely determined by the initial data ([33]), and 
moreover their spatial support will not increase with time ( [7] 133] [56] ). 

Even in the case when one resorts to continuous initial data which are strictly positive throughout 
n, in which in fact unique classical solutions exist for any p > 0, the value p = 1 corresponds to a 
critical strength of degeneracy in that for p < 1, after an appropriate waiting time, all solutions will 
enter the cone := {(p ; —)• M | (p(x) > cdist(x,9n) for all x G and some c > 0} ([6]), thus 

reflecting a diffusion-driven effect generalizing the Hopf boundary point property in non-degenerate 
diffusion processes. In the case p > 1, however, solutions to (|I.2p emanating from initial data which 
are suitably small near dQ will never enter l|55jl. 

Now in (|1.II) . this degenerate diffusion process interacts with a spatially nonlocal source which is such 
that unlike in large bodies of the literature on related nonlocal parabolic equations (|lQ]), even basic 
questions concerning local solvability appear to be far from trivial: Indeed, in light of an expected loss 
of appropriate solution regularity due to strongly degenerate diffusion, even for smooth initial data it 
seems a priori unclear whether solutions can be constructed which allow for a meaningful definition of 
the Dirichlet integral for positive times. This is in stark contrast to most nonlocal parabolic 

problems previously studied, in which either diffusion is non-degenerate and hence such first-order 
expressions are controllable by bounds for solutions at least for small times, such as e.g. in the 
semilinear problem 


ut = Au + u^ |Vnpdx 

studied for m > 1, r > 0 in m). or the nonlocal terms involve only zero-order expressions which 
thus in a natural manner also in cases of degeneracies as in (jl.2p allow for local theories based on 
extensibility criteria in L°°(D) only (see [101144] and also the book [40] ). 

Main results. Previous mathematical studies on the PDE in (|l.ip have concentrated on analyzing 
self-similar solutions only. In [23], the authors constructed self-similar solutions in the case D = M, and 
in [38] the same could be achieved in the multi-dimensional case D = with N > 2. More recently, 
the authors in |39j investigated the existence of self-similar solutions in the one-dimensional case in 
a closely related problem in which the Laplacian is perturbed by a time-dependent term containing 
the first derivative as well; all these self-similar solutions are shown to be singular and to approach 
Dirac-type distributions as t \ 0. 

The goals of the present work consist in developing a fundamental theory of local solvability for 
and in providing a first step toward an understanding of the qualitative solution behaviour. In 
order to formulate our results, let us concretize the specific setting within which (II.Hi will be studied 
by requiring that throughout the sequel, D denotes a bounded domain in > 1, with smooth 

boundary, and by introducing the solution concept that we shall pursue as follows. 

Definition 1.1. Let T € (0,oo]. By a weak solution of il.l]) in LI x (0,T) we mean a nonnegative 
function 



u G X [o,r)) n LU%T);W^’\Q)) 


with Ut G LL(^x[0,T)) 
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which satisfies 

rT 


— f fu(ptdxdt+ f fvu ■V{uip) dxdt = fuo(p{-,0) dx + f ( fuipdx') ■ ( f\V ul"^ dx') dt (1.3) 
JoJn JoJn Jn Jo^Jn ^ ^Jn ^ 

for all^e C^{n X [0,T)). 

A weak solution u of in Q x (0,T) will be called locally positive if ^ ^ [Oj^])- 

Remark 1.2. Sinceu G T); lkQ^’^(r2)) andut G L‘i^fiQx[0,T)) imply thatu G C^{[0,T)-, L‘^{Q)), 

111.A) is equivalent to requiring that n(-,0) = uq, and that 

f [ utipdxdt+ f f Vu ■ V{u(p) dx dt = f ( f u(pdx\ ■ ( f iVul^dx^dt (1-4) 

Jo Jn Jo Jn Jo ^Jn '' ^Jn ' 

holds for any ip G x (0, T)). 

In order to construct such locally positive weak solutions, we shall assume that the initial data satisfy 


(HI) no G L°°{n) n IVo’^(O) and 
(H2) no > 0 and ^ G as well as 

(H3) there exists L > 0 such that ||no||$ < L. 

Here and below, for a measurable function v. H —>• M we have set 


Pll<i-,oo ■= ess sup 


where G C^(H) denotes the solution to 


— A<b = 1 in 


4*130 = 0. 


(1.5) 


Note that according to the Hopf boundary point lemma, requiring ||no||,j, to be finite is an equivalent 
way to ask for the possibility of estimating no by a multiple of the function measuring the distance of 
a point to dkl. 

In this framework, the first of our main results indeed asserts local existence of locally positive weak 
solutions, along with a favorable extensibility criterion only involving the norm of the solution in 
L°°(H). 

Theorem 1.3. Let uq satisfy (H1)-(H3). Then there exist T^ax G (0,oo] and a locally positive weak 
solution u to ([HD in Q X {0,Tmax) which satisfies 

either Tynax = oo or limsup ||n(-, t)||^oo(Q) = oo, (1-6) 

t/'Tmax 

and which is such that for each smoothly bounded subdomain fl' CC H there exists Cq' > 0 with 


[ |Vn(-,t)p < [ |VnoP • exp 
Jn Jn 


1 


sup / n(-,r) 


2Co/ \T&(o,t)Jn 


(/)lnn(-,t)— / (/>lnno + 


'n' 


in' 


0 Jn' 


as well as 


|«(-,i)ll<i.,oo ^ ™ax < ||no||$^oo > sup \Vu{x,t)\'^ dx 
I Te(o,t) J n 


(1.7) 


( 1 . 8 ) 


foT Q,.€. t E if^iTYnax) • 
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We emphasize that the extensibility criterion ()1.6|) particularly excludes any gradient blow-up phenom¬ 
enon in the sense of finite-time blow-up of Vn despite boundedness of u itself. Indeed, the occurrence 
of unbounded gradients of bounded solutions appears to be a characteristic qualitative implication of 
various types of interplay between diffusion, possibly degenerate, and gradient-dependent nonlineari¬ 
ties ([aillEIlllg]). 

A natural next topic appears to consist in deriving conditions on the initial data which ensure that 
the solutions found above either exist for all times, or blow up in finite time. Here in view of the 
essentially cubic character of the production term in (|l.ip it is not surprising that this may dominate 
the smoothing effect of the merely quadratic-type diffusion term when the initial data are suitably 
large in an adequate sense; precedent works indicate that indeed such intuitive considerations are 
appropriate in related non-degenerate and degenerate parabolic equations with local reaction terms 

([MIllllSlEH] ). 

As a remarkable feature of the precise structure of this interplay in (11.11) . we shall see that actually 
a complete classification of all initial data in this respect is possible, exclusively involving the size 
of the total initial mass m := JqUq as the decisive quantity: In fact, the second of our main results 
identifies the value m = 1 to be critical with regard to global solvability, and moreover gives some 
basic information on the asymptotic behaviour of solutions. 


Theorem 1.4. Let uq satisfy (H1)-(H3), and let u and T^ax denote the corresponding locally positive 
weak solution of dni), as well as its maximal time of existence, provided by Theorem 11.31 
(i) If f^uo < 1, then T^ax = oo and 


/ u{x, t)dx ^ 0 as t ^ oo. 

Jn 

(a) Suppose that JqUq = 1. Then T^ax = oo and 

/ u{x, t) dx = 1 for all t > 0. 

Jn 

(in) In the case f^uodx > 1, we have T^ax < oo and 


lim sup 

i/^Tmax 


/ u{x,t) dx = oo. 

Jn 


Remark 1.5. The statement (ii) of Theorem \l.4\ says that if the initial data uq is a probability measure 
then we have conservation of probability in time. This is actually a desired feature of the replicator 
dynamics model described by (ini), since u{-,t) stands for a probability distribution of the state of 
some population of players, see also Appendix A. 


In the situation of Theorem 11.41 (hi) when finite-time blow-up occurs, understanding the solution 
behaviour near the respective blow-up time necessarily requires to describe the set of all points where 
the solution becomes unbounded. Accordingly, we shall next be concerned with the blow-up set 


= S X € H 


there exists a sequence C H x (0, T^ax) such that 

Xk x,tk ^ Tmax and u{xk, —)• oo as A: —)• oo > 


of exploding solutions. In numerous related equations, involving either linear or degenerate diffusion, 
blow-up driven by local superlinear production terms is known to occur in thin spatial sets only, in 
radial settings typically reducing to single points ([a ttzi [12]), with only few exceptional situations 
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detected in the literature which lead to regional or even global blow-up, thus referring to cases in 
which |,^| > 0 or even = Q. (cf. [lU [THl [23 SHI IM] , for instance). In cases of sources which at least 
partially consist of nonlocal terms, blow-up in sets of positive measure may occur if the relative size of 
a possibly contained local contribution at large densities is predominant, as compared to the strength 
of the respective diffusion term ([niisoi i 32 i sa sa EH). 

Our main result in this direction will reveal that any of our non-global solutions in fact blow up globally 
in space, thus indicating a certain balance in the competition of diffusion and nonlocal production in 

O)- 

Theorem 1.6. Suppose that f^uodx > 1, and let u denote the locally positive weak solution of (11.111 
from Theorem \1.3l Then u blows up globally in the sense that its blow-up set satisfies = Q. 

The outline of the paper is as follows. In Section [2] we introduce an approximate sequence of non¬ 
degenerate problems and derive some estimates for their solutions u^. Here one key step toward the 
existence proof will consist in deriving the associated approximate variant of (12.351) fLemma l2.6l) . wich 
will rely on an energy type argument combined with an analysis of the functional InUs{-,t) for 
n' CC n, f > 0 and appropriate fi. Another important observation, based on an integral estimate 
involving certain singular weights (cf. Lemma 12.51 and in particular (I2.3ip l. will reveal that the 
functions Vwe enjoy a favorable strong compactness property with respect to spatio-temporal L^- 
norms (cf. (j2.44p ). rather than merely the respective weak precompactness feature obtained from 
corresponding boundedness results. In Section [3] we study an ODE problem associated with the 
evolution of the total mass of the solution, and in dependence on whether this total mass initially is 
equal, less or greater than 1, we prove global existence and conservation of the total mass, convergence 
to zero total mass and finite-time blow-up, respectively. Finally, in Section 0] we concentrate on the 
latter case and examine the corresponding blow-up set of the solution, and we actually prove that any 
such blow-up occurs globally in space. 


2. Weak solutions: existence and approximation 


Following an approach well-established in the context of degenerate parabolic equations, we aim at 
constructing a solution u to m as the limit of solutions to certain regularized problems. For 
this purpose, let us fix a sequence C (0,1) such that Sj \ 0 as j oo, and a sequence 

{uoe)e=ej C with the properties 


uoe > e in D, = e on 5D, Aiioe = — |VrtoeP on 5D for all e G 

Jn 

and 

limsup lluoe - e||<i,,oo ^ 

£=£j\0 

with L > max I J^lVttoP, lliioll^ oo}> (^3)) 

for any compact set iL C D there is Ck > 0 such that liminfinfMo£ ^ Ck, 

£\0 K 


and such that moreover 
and 


mo£ —>■ in VF^’^(D) as e = \ 0 



for all £ G (ej)jgN. 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 
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A necessary first observation is that such an approximation actually is possible. 

Lemma 2.1. Let uq satisfy (H1)-(H3). Then there is a sequence ^ having the 

properties 

Proof. By modification of the usual mollification procedure (cf. [581 Section 13]) commonly employed 
to obtain ()2.4I) it is possible to obtain the other properties as well. More precisely, we set 

U 0 e = e + C{1 - p)4> + p{tp + ai?), 

where p € C'^(n) is a mollified version of uq (after “locally shifting uq towards the interior of the 
domain”), p G C^(fi), 0 < p < 1, such that the supports of Vp and ip are disjoint, 0 < i? G with 
J^'d = 1 (in order to adjust (I2.5p i. is the solution to —= 1 in fl, = 0 on dLl (for achieving the 
third property in (|2.ip 1. and C and a are appropriately adjusted constants, depending on e as well as 
several different integrals containing the functions $, p, i?, their gradients, and uq- 
For a slightly more detailed version of the proof, we refer the reader to the appendix. □ 


For e G (ej)jgN, we consider the regularized problem 

Uet = + Ug • Jj^lVUg] ), 

< Ue{x,t) = e, 

. Ue{x,Qi) = UQe{x), 


X G 12, t > 0, 

X G dLl, t > 0, 
X G 12, 


where 


pJ.z) 


:= min 



for ^ > 0. 


( 2 . 6 ) 


Lemma 2.2. For all sufficiently small e G (ej)jgN; problem \2. h|) has a unique classical global-in-time 
solution Us G C^’^(12 x [0,oo)). 


Proof. To prove the uniqueness statement for all e, we assume that both ui and U 2 are classical 
solutions of (12.61) from the indicated class in 12 x (0, T) for some T > 0. Then w := ui — U 2 satishes 
re = 0 on dQ and at 2 = 0, and 


Wt = UiAw + Att 2 ■ W -\- Ps\ I |Vtt 2 r ) ■ W -\- Ui 


Pe\ 


\Vuif)-pJ / IV 1 X 2 I 


(2.7) 


for t G (0,r). Now given T' G (0,T), we can find a constant M > 0 such that ui,|Vui|,tt 2 and 
IVU 2 I are bounded above by M in 12 x (0,r'), since ui,U 2 are classical solutions. Thus, by Holder’s 

inequality and the pointwise estimate |Vui| — |Vrt 2 | < |V(tti — U 2 )|, we obtain 


Ps[ljVui\‘^)-ps{jjVu2\‘ 


<\\Pe\ 


L°°((0,oo)) 


(|Vni|2-|Vn2|2) 


< 


f |Vtti|-|Vu2| • (|Vtti| + IVU 2 I) 
Jn 

<2M / |Vu;| 

Jn 

<2M|12|I • ( [ iVixp''" 


( 2 . 8 ) 
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for all t G (0, r'), because ||pe ||loo((o^oo)) < 1- Upon multiplying (j2.7l) by w and integrating over Q we 
see that for t G (0, T') 


1 d r 

2 dt Jq 


= 


(2.9) 


< - 


/ Ui^WW + / /^U 2 + / Pei / |Vu 2 |' 

Jo, Jfi Jq \Ju 

/ ui\Vw\' 

Jn 

+ [w^Pei [ |Vu2p) + [ |u;|ui pj f |Vuip^ - Ps( [\VU 2 \‘ 

JQ \Jq J Jq \Jq J \Jn 

Together with Young’s inequality, (12.8p and the facts that ui > e and pj^s) < ^ for all s > 0, this 
entails 

1 d 


p 4 j - p,^JJVu2\- 


— / VuiVuiw — 2 / wVwVu2 
Jq Jq. 


2 dt 


I 

Jq 


< -e 


[\^w\^ + 7 [ |Vu;|2 + i /'u;2|Vui|2 + ^ [ |Vu;|2 + - [ w^\Vu2(^ 
Jq ^ Jq ^ Jq J Jq £ Jq 

1 

+ - [ w^ + 2M\n\'2 ( [ |Vrc|^ 

£ Jq \Jq 


\W\Ui 


for t G (0,T'). The choice of M now ensures that 


1 d 


JVP 


w^<-^ / |Vu;|^ + — w^ + 


8M2 


2dt 


£ Jq 


£ Jq 


+ 


i/ 

£ Jq 


w 


+ 2M|f2|2 ( / |Vu;|' 

Iq 


/ / ^^1 

'Q Jq 


< - 


[ iVinl 

Jq 


, 9M2 -h 1 

^ + 


4./o' e 

for t G (0,T'), so that (12.101) finally turns into 


2 ^Jq £ 


/ 

Jq 


w 


( 2 . 10 ) 


2 dt 


I 

Jq 


w < 


Id f 2 /9M^ + 1 


+ 


w 


e e 

for all t G (0,r'). 

Integrating this ODI yields that rc = 0 in fl x (0, T') and hence also in 12 x (0, T), because T' <T was 
arbitrary. 

It remains to be shown that for all T > 0, (12.6j) is classically solvable in 12 x (0, T) provided e is 
sufficiently small. To this end, fix T > 0 and let e G (ej)jeN be so small that /j^|Vrt£p < which 
is possible due to (12.41) . By [281 Thm. V.1.1], there are Ki > 0 and 0 > 0 such that any classical 
solution w to the problem 

wt = wAw + c{x,t) in X [0,T], 'w\qq = s, w{-,0)=uo£ 
with c G L°°{Q X (0,T)) fnlfilling 0 < c < ^ ||7i0£|lioo(n) which in addition obeys the estimate 

T 

£ <w < ||uoe||oo eT satisfies 

ll'R’II a e — < Kl- 

" "C®’^(f2x[0,T]) “ 


( 2 . 11 ) 
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Fix (5 > 0. Corresponding to 9, Ki and 6, there is K 2 such that any solution w to 

wt = a(x, + h{x, t) in O x [0, T], w\dn = c, 0) = uoe 

for some a G x [0,T]) having the properties a{x,t) = e for {x,t) G 90 x [0,T], e < a < 

T 

lluoell^oc eT, < ^1 and continuous b with b{x,0) = bo £ R, ||6| 

application of m Thm. 7.4] to tc — uoe — tbo fulfils 


< by an 


R’ — < K2. 

' "c^+'*’^(r2x[0,T]) “ 

s 


( 2 . 12 ) 


With this in mind, in the space X = x [0,T]) we consider the set 


5 := u G X 


u > e in O X (0,T),u(-,0) = uoe and ||u|| ,,, s — < 


K 2 }, 


which is evidently closed, bounded, convex, and compact in X. For each v £ S, the definition of pe 
implies that 


/(t):=p,(^|Vu(-,t)|2), tG[0,r], 


(2.13) 


defines a nonnegative ^-Holder continuous function / on [0,r]. The choices of f,S and £ show that 
/(O) = /f^jVttoeP and thus (12.ip ensures that the compatibility condition of first order is satisfied. 
Therefore, the quasilinear, actually non-degenerate parabolic problem 


Uet = UeAUe + f{t)Ue, 
Us{x,t) = £, 


X £ ^}, t > 0, 

X £ dQ, t > 0, (2.14) 

. Us{x,0) = UQe{x), X £Q, 

possesses a classical solution Ue £ C'^’^(n x [0,r]) by [28l Thm V.6.1], which, by comparison, satisfies 

T 

e < tie < I|ji 0 £||i,<-( 0 ) • e“ inflx(0,T), (2.15) 

because u{x,t) := e and u{x,t) := ||uo£||l°°(o) • are easily seen to define a sub- and a supersolution 
of (|2.14p . respectively. 

We now introduce a mapping F: S ^ X hy setting Fv := u^, where solves (I2.14h with (|2.13p . 

Then defining c{x,t) := Us{x,t)f{t), x £ Q,t £ [0,T], this function satisfies ||c||gQ < ^ ||uoe|loo®^ 
accordingly, as stated in (|2.11l) above, jlFujl g e < Ki for any v £ S. 

Using a{x,t) := {Fv){x,t) and b{x,t) := {Fv){x,t) ■ f{t), we see that, again, the above considerations 
are applicable and s— < K 2 for any u G S' by (I2.12F In particular, we observe that 

FSdS. 

Furthermore invoking [28] IV.5.2], we can conclude the existence of /c > 0 and > 0 such that 


IlFujl 


C^+'*’i+^ (Dx[ 0 ,T]) 


< k + II'“0£|Ic2+^(dx[0,t]) + ^ -^3 (2.16) 

for all V £ S. To see that F is continuous, we suppose that {vk)k£N C S and v £ S are such that 
Vk^ V in X. Then fk{t) := /^|V?;fc(-, t)p) satisfies 

fk^f mC°([0,r]) 


(2.17) 
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as k ^ oo, with / as given by (j2.13jl . By (j2.16p and the theorem of Arzela-Ascoli, {Fvk)k£^ is 
relatively compact in x [0,T]), and if fej —>• oo is any sequence such that Uki ■= Fvki converges 

in X [0, r]) to some w as i ^ oo, then in 

dtuki = ukiAuki + fki{t)uki, X e n, t G (o,r), 

we may let fcj —)• oo and use (|2.17l) to obtain that w is a classical solution of (|2.14p . Since classical 
solutions of (I2.14P are unique due to the comparison principle, we must have w = Fv. We thereby 
derive that the whole sequence {Fvk)keN converges to Fv and hence conclude that F is continuous. 
Therefore the Schauder hxed point theorem asserts the existence of at least one Ug € S for which 
Us = Fus holds. Since such a fixed point obviously solves (12.61) . the proof is complete. □ 


The basis of both our existence proof and our boundedness result is formed by the next two lemmata 
which provide useful a priori estimates for in terms of certain presupposed bounds. The first lemma 
essentially derives a uniform pointwise bound for from a space-time integral estimate for |Vu£p. 


Lemma 2.3. For all M > 0 and B > 0 there exists C{M,B) > 0 with the following property: If 

rT 

(2.18) 


uoe < M in It and f f |Vu£p < B 

Jo Jn 


holds for some e G and T G (0, oo] then we have 

Ue < C{M,B) in H X [0,T). (2-19) 

Proof. Our plan is to use a separated function of the form 

u{x, t) := y{t) ■ (M -|- 4>(x)), x G O, t G [0, T), (2.20) 

as a comparison function, where M is as in the hypothesis of the lemma, G (7^(12) is the solution of 
and y denotes the solution of 


y' = + {f{t) +1) • y, (o, r), ?/(o) = i, 

with f{t) := /^|Vu£(-, t)p. In fact, it follows from (I2.2ip that z := i is a solution of z' 
2 :( 0 ) = 1, and hence given by 

Z{t) = e- Jo ^ f' g- /J f(a)da-it-s) ^ ^ ^ 

Jo 

We claim that 

1 < y{t) < 6'®+^ for all t G (0, r). 

To see this, we note that if t G (0,T) satisfies t < 1, then (I2.18P implies 

z{t) > e-Jo > e-^-* > e-^-\ 

whereas if t G [1,^) then again (I2.18P shows 

z{t) > f e-JsJ^^^'^^-J-^Us> f e-^-J-^Us 

Jt-i Jt-i 

> f e-^-^ds = e-^-\ 

Jt-i 


( 2 . 21 ) 


( 2 . 22 ) 
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This yields the right inequality in (I2.22|) . while the left immediately results from an ODE comparison 
of y with y{t) = 1, because y' + y'^ — {f{t) + l)y = —f{t) < 0. Consequently, since > 0 in D, the 
function u defined by (j2.2U|) satisfies 

u{x, 0) = M + ^(x) > M > Ue{x, 0) for all x G D 

due to (|2.18l) . and on the lateral boundary we have 

u{x,t) = y{t) ■ M > M >e for all x G 5D and t G (0, T). 

Moreover, 

ut - uAu - f{t) -u = y' ■ {M + ^) + y^ ■ {M + ^) - f{t)y • (M + $) 

= y • (M + ^>) 

> 0 for all a; G D and t G (0,T), 

whence the comparison principle ensures that <u in Qx (0,T). In view of (|2.22l) . this entails that 

Ue{x,t) < ■ (M + 11 $11^00(^7)) for all x G 0 and t G (0,f), 

so that (I2.19P is valid upon an obvious choice of C = C{M, B). □ 

Next, the fact that solutions of (12.6p cannot blow up immediately can be turned into a quantitative 
local-in-time boundedness estimate in terms of the norm of the initial data in L°°(D) n 
Moreover, our technique at the same time yields an estimate involving integrals of Uet and Vtte, as 
long as Ue is appropriately bounded. 


Lemma 2.4. i) For all M > 0 there exist Ti{M) > 0 and C\{M) > 0 such that if 


uoe < M in Ft and / |VttoeP < M (2.23) 

Jn 

hold for some e G then 

Ue < Ci{M) in Fix [0,Ti(M)). (2.24) 

ii) For each M > 0 and T > 0 there exist T 2 {M) G (0,T] and C 2 {M) > 0 such that whenever 
£ G (ej)jgN is such that 


Ue < M in Ft X (0,T) and / |VrtoeP < M 

Jn 


are satisfied, then 


rT2{M) r,..z i- 

/ /^+ sup \VUei;t)f <C 2 {M). 

Jo JnUe te{o,T2{M)) Jn 


(2.25) 

(2.26) 


Proof, i) We multiply (|2.6I) by ^ and integrate by parts, use that Uet = 0 on dQ, and apply Holder’s 
together with Young’s inequality to see that 



(2.27) 












A NON-LOCAL DEGENERATE PARABOLIC PROBLEM 


11 


for all t > 0, because Pe{z) < z for all z > 0. Hence, 



Using the Poincare inequality, we obtain 


(2.28) 


<ci- +lj 

with a positive constant ci independent of e G (ej)jgH £ (0,1) and t > 0. Therefore, (12.281) yields 

+ (/JV„.p)^ (2.29) 

which in particular implies that y{t) := J^|Vu£(-,t)p satishes 

y'{t) < c{^/y + l)y^ for all t > 0. 


Hence, if we let z denote the local-in-time solution of 

z'{t) = c{-sfz -t- 1 ) 2 ^, 
z(0) = M, 


f > 0, 


with maximal existence time > 0, then due to (j2.23p and an ODE comparison we have y < z in 
(0, Tx). Defining Ti(M) := for instance, we obtain from this that < z{Ti{M)) for 

all t G [0, ri(M)), whereupon (|2.24p now results from Lemma 12.31 

ii) If the first inequality in (I2.25P holds then (I2.28P entails that y as defined above even satisfies the 
nonlinear ODI 


y'{t) < M\Q\y‘^ for all t > 0, 

whence we have ^ ^ with T 2 := min{r, 1/(M^|D|)}, by 

the second inequality in (j2.25p . Inserting this into (I2.29p again and integrating over (0,T2) proves 

UM- □ 


When constructing the solution u of dEU) as the limit of solutions of (|2.6p . it will be comparatively 
easy to obtain the approximation property Vug —)• Vu in the sense of x [0, r))-convergence. 

For handling the nonlocal term in the equation, however, it seems appropriate to make sure that also 
^|Vu,|2^^|Vu|2inLU[0,r)). 

In order to achieve the latter we exclude certain boundary concentration phenomena of in the 
following sense. 


Lemma 2.5. For any T>0, C>0, M>0 and 5 > 0, there is K = K{M,C,T,6) CC D and rj > 0 
such that whenever e G (ej)jgN is such that e < y and 


we have 


sup / |Vue(f)p < C and < M, 
te[o,r] JQ. 


r [ ivue 

Jo Jn\K 


< 5 . 


(2.30) 
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Proof. For q G (0,1), we multiply (j2.6p by nf ^ and integrate by parts to obtain 


1^ 
q dt 


[ ul= [ ujduUe- f qul ^\VUe\‘^+ f ul Pei f \VUe 

Jn Jan Jn Jn \Jn 


where we can use dyUe < 0 on 90 and integrate with respect to time to derive 


- r [ur^\vu^ 

Jo Jn 


el < — 


1 

Q Jn 


ul{T) + 


1 

<1 Jn 




\Vue\^ =: C(r) 


(2.31) 


for all e > 0 satisfying (I2.30D . which gives control on |Vttep whereever Ug is small - which is the 
case near the boundary, as we ensure next: In order to lay the groundwork for the corresponding 
comparison argument, note that by (I2.3njl . 


Uet = UeAUe + UePs 


[ NUe 

Jn 


^ 'a.eAUe “1“ CUe^ — C, Ue(0) — noe* 


(2.32) 


Fix p > 0 such that ^C(T) ^ ^ solve (II.5F Choose A> C such that + p > uqs for all 

0 < £ < p, which is possible due to condition (j2.2j) . Then u := A^ + p satisfies 

ut = 0 > — (A<I> + p)A + (^$ + p)C = uAA^ + Cu = uAu + Cu. 

As long as £ < p, also u|go > tie bo holds and furthermore 

m( 0) > rtoe- 

Therefore, by the comparison principle, we obtain u > Ue- 
Now choose iF CC fl in such a way that 

A^ < p in 12 \ iF. 

This entails Ue <u = A^ + p < 2p in D. \ K. Then 

rT 


f [ |Vn,p = [ [ ut-^\Vue\\l-'^ 

Jo Jn\K Jo Jn\K 

< (2r?)i-'^ r [ 

Jo Jn\K 

< {2p)^-^ r [ 

Jo Jn 


e\ ^ 


{2p)^-<^C{T) 


q 


□ 


by virtue of (I2.3ip . 

We are now ready to prove that the Ue in fact approach a weak solution of (II. ip that is locally positive 
in the sense of Definition 0 Before we do so, however, we prepare the following estimate for Ue that 
will be useful in proving assertions about the blow-up behaviour of u. 

Lemma 2.6. Let 12' CC LI be a domain with smooth boundary. Assume also that (p denotes the 
solution to —A(f> = 1 in 12', 4>\qq' = 0. Then there exists Cq/ > 0 such that for each e G (ej)jgH and 
any t > 0 the solution Ue of (12.6p satisfies 


\^Uei;t)\^ < 


iVnoepexp 




sup 
, Te(o 


3 / Ueir) 

,t)Jn 


/ 4>lnuei-,t)- / blntto£+ / / Ue 

In' Jn' Jo Jn' 



(2.33) 














A NON-LOCAL DEGENERATE PARABOLIC PROBLEM 


13 


Proof. As Ust = 0 on dQ, similarly to (|2.27p . multiplying (|2.6p by ^ and integrating over 0 yields 

f = [ Uet^Ue+ f UetPei [ 

Jn u jfi Jq \Jq j 

IVRel^A / UetPei / IVlig 

Jq \Jq 


1 d r 

2 dt Jq 

After rearranging, by Holder’s and Young’s inequalities and the definition of p^ this entails 




.1 < -2 


f 

Jn Ue 


+ 2 


Uet 


n\V^‘ 


1 

2 \ 2 


[ 

Jn 


Ps 


In 


IVUe 


< - 


< 


2 f ^ + 2 [ ^ + 1 [ UePe( [ |Vn, 

Jn Jn 'H'e 2 Jq \jQ 

^ j f on (0,Oo). 


In J Jn 

This looks like a quadratic differential inequality for y{t) = and at first does not seem helpful 

for obtaining an estimate for this quantity. Therefore we shall split the respective quadratic term and 
apply Gronwall’s lemma to y'{t) < g{t)y{t), where 


which leads to 


y{t) < y{0) exp f g{T)dT for all t > 0. 

Jo 


(2.34) 


In this situation, however, we are left with a term fg in the exponent and we prepare an 

estimate for this in the following way: With J) as specified in the hypothesis, we let Cq' = 4> > 0. 

Multiplication of ()2.6I) by ^ and integrating over Q' then gives 


— / Inuej) = / 

m Jqi Jq/ 


AUe4> + 


/ J^Pei [ IVUe 
Jn' \Jn 

/ UeAcj)- / d^Ue4>+ / UedyCf + Cq/P el 

Jn' Jdn' Jan' \j 


Wup 


on (0, oo). 


Taking into account the definition of f and its consequence < 0 = J^qq'., we infer that 


^ [ (j)lnUe>- [ Ue + CQ/pe 

Jn' Jn' 


iVriel'’' ) on (0, oo). 


Therefore 


Pe 


^ cj. 



0 Jn' 


Us+ (t>lnup{t) - / 4>lnuoe 

Jn' Jn' 


for any t > 0, and we can conclude from (j2.34l) that 


f |V'U£;(t)p < [ iVRoel^exp sup f u{t) ( f f Up+ [ (plnup{t) - [ cflnuoe 

In Jn ‘^'-'n' Te(o,t) Jn \Jo Jn' Jn' Jn' 


for all t > 0. 


□ 
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Another useful piece of information is that a condition like (H3) remains satished for any t > 0. 


Lemma 2.7. Let T > 0, M > 0 and e G (ej)jgN be such that ||uoe — e||$ < oo. Then any solution 

Ue of (1^ which satisfies 


already fulfils 


|Vu,(t)|2 < M 


for any t G [0, T] 


U. - £ 


I'LjCxd — 


maxjM, ||uo£ 



Proof. Let C = max{M, ||Moe ~ £||$ ooi consider u := + e with <i> as in (jl.5p . Then ut = 0 > 

(M —C)(C'<h + e) = uAu + Mu, whereas u^t = "UeAue+ rt£/Oe(/f^|Vttep) < UeAus + Mu^. Additionally 
u\dn = £ = Us\dn and u{x,0) — £ = C^{x) > <h(x) ||ttoe — e||$oo A uofix) — £ and therefore the 
comparison principle [S2] asserts that <u and hence implies the claim. □ 


With this information at hand, we can proceed to the proof of convergence of the Ue to a solution of 
(jl.ll) that still satisfies an inequality like (j2.33p . 


Lemma 2.8. Suppose that uq satisfies (H1)-(H3). Then there exists T > 0 depending on bounds on 
II'“o|Il°°(d) II^'^o|Il 2 (q) and a locally positive weak solution u of in Ll x (0,T). This solution 

can be obtained as the a.e. pointwise limit of a subsequence of the solutions of h2. fij) as £ = £j \ 0, 
and for any smoothly bounded subdomain O' CC O there is Cq/ > 0 such that 


[ |Vu(-,t)l2 

Jn 

< f iVuoPexp 

Jn 

as well as 


for a.e. t £ (0,T). 



(2.35) 


(2.36) 


Proof. We set Mi := max{||uo||Loo(Q) + 1, /q|VuoP + l| and let Ti = Ti(Mi) and ci = C'i(Mi) be as 
in Lemma 12.41 ih Then this lemma states that Ug < ci in 0 x (0,Ti) for all £ G (ej)jgN- Accordingly, 
corresponding to M 2 = max{ci, J(^|VuoP + 1}, Lemma [2.41 hi provides T = T 2 (M 2 ) G (0,Ti) and 
C 2 = ( 1 ^ 2 (M 2 ) > 0 such that 

[ [ — + sup / |VUe(-,t)|^ < C 2 (2.37) 

Jo Jn'^e t£{o.T)Jn 

for all £ G (ej)jgN, which by < ci can be turned into a uniform bound on ||wt|lL2(Qx(o T)’ from 
which it follows by means of the fundamental theorem of calculus that after possibly enlarging C 2 , we 
also have 


\u. 


'C2([0,T];L2(D)) 


< C2 


(2.38) 


for such e. 

In order to prove a uniform estimate for from below, locally in space, we follow a standard compar¬ 
ison procedure: we pick any smoothly bounded domain 14' CC 12 and let (f) G C‘^{Pt') solve —Acf = 1 
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in ri' with = 0. Then the lower estimates in ()2.3|) guarantee that 

uoe{x) > C 3 (n')(/)(x) for all x G n' (2.39) 

holds with some appropriately small C 3 (n') > 0. Letting y{t) := > 0, denote the solution 

of y' = —with y( 0 ) = 03 ( 0 '), we thus find that u{x,t) := y{t)(f){x) satisfies u<u^ on the parabolic 
boundary of 12' x (0,oo). Since 

Uf- — uAu = y'(j) + y'^cf) = 0 in 12' x (0, 00 ) 


and 

Uet - UeAUe = ■ Pe(^ j | >0 in 12 X (0, oo), 

we conclude from the comparison principle (see [52] for an adequate version) that u<Ue and thus, in 
particular, that for each compact K <Z Vt and T' > 0 there exists a suitably small C 4 {K,T') > 0 such 
that 

Ue > ca{K, T') in K X (0, T') (2.40) 

holds for all e G Now the estimate < ci, (I2.37|) . (12.381) and (I2.4n|) along with standard 

compactness arguments allow us to extract a subsequence of (ej)jgH and a function u : 

12 X [0, T] —)• M such that 


Ue ^ u 

in C^([0, T); L^(12)) and a.e. in 12 x (0,T), 

(2.41) 

Vti£ — 

^ Vu 

in Tfo^(i2 X [0,r)) and 

(2.42) 

Uet - 

Ut 

in L^{VL X (0,r)) 

(2.43) 


as e = £jf, \ 0. From (I2.4ip . the inequality Ue < Ci and (I2.4U|) . we know that u < ci a.e. in 12 x (0, T) 
and u > C 4 ^{K,T) a.e. in iL x (0,T) whenever K CC 12. Moreover, since tte — e vanishes on 512, (|2.42p 
implies that u G L^((0,T); (12)), so that u fulfills all regularity and positivity properties required 

for a locally positive weak solution in 12 x (0,T) in the sense of Definition ll.il 

In order to verify that rt is a weak solution of dni) it thus remains to check (jl.4l) . To prepare this, we 
claim that in addition to (j2.42p . we also have the strong convergence properties 

Vug —>• Vu in ^ in f2 x (0,T) (2.44) 


as well as 


|Vue(x, •)pdx — )■ f \Vu{x, dx inL^((0, r)) 

i Jn 


(2.45) 


as e = £ji^ \ 0. To see (I2.44p . we let iF CC D be given and fix a nonnegative 4) G (^(^^(ri) such that 
i/; = 1 in iL. Then 

rT r rT 


[ [ IVtte-Vup 

JQ J K 


< 


/o Jn 

rT 


/o JQ. 


JQ 

/ V(u£ — u) ■ Vui; ■ 'Ip + / Vu • V(u£ — u) ■ Ip 

JQ Jo Jq 


=: /i(e) + / 2 (e) for all e G (ej)jeN, 


(2.46) 


where / 2 (e) —)■ 0 as e = \ 0 by (12.421) . Using the equation for Ue, however, after an integration by 

parts we find that 


h(e) 




I 

JQ 


{Us — u)VUe ■ V'lp 
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T T 

= -f — [(Ue-u) ■ Pe( [ -'ll; 

Jo Jn '^e Jo Jn ^Jn ^ 

- f [{u£- u)Vus ■ Vi/j 

Jo Jn 

=■ hii^) + for all e G (ej)jgN. 


Due to (j2.4ip and (j2.42p . we have / 13 (e) —)• 0, and (j2.4ip together with (I2.37p and Holder’s inequality 
imply that 


I^i 2 (e)| < 



• IIV’IIl2(o) 0 


as e = Ej^ \ 0, where we again have used the fact that pj^z) < z for any z > 0 and all e G (ej)jgN. 
We now use Holder’s inequality and the local lower estimate ()2.40p . which in conjunction with (I2.37P 
yields 


|Ai(e)| < 
< 



IIV’IIl°°(o) 

(C4(suppi/;,r))2 




^ 0 


as e = Ej^. \ 0, by (I2.4ip . Altogether, we obtain that /i(e) —)• 0 and hence, by (I2.46p . that Vug —Vu 
in L^{K X (0,T)) as e = e^j. \ 0 for arbitrary K CC D. 

Having thus proved (j2.44p . with the aid of Lemma 12.51 we obtain (I2.45p as a straightforward conse¬ 
quence: 

Given (5 > 0, we let K = K{ci,C 2 ,T, |) and p > 0 he the set and the constant provided by Lemma 
12.51 and employ the convergence asserted by ()2.42p to choose ko G 'N such that for all k, I > ko we 
have ^ Then for all k,l > ko, 




< r [ ||Vu,j2-|Vu,j2|+ r f 

Jo Jk Jo Jn\K 

+ r j \vusf 

Jo Jn\K 

6 6 6 

< K + T + T 


and thanks to the completeness of L^((0, T)) we obtain (j2.45p . We can now proceed to verify that 
(11.41) holds for all if G X (0,T)). To this end, we multiply (|2.6p by y? G x (0, T)) and 

integrate to obtain 



Uet(P + 




UeVtte • Vip 


In 


Here, as e = \ 0 we have 



• 7^- 
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by (I2.4.S|) . whereas (I2.44|) and (j2.4ip allow us to conclude that 


r [ r [ |Vn|‘ 

Jo Jq Jo Jn 




and 


/ / UeVUe'Vy?—)■ / / uS/u-Vip, 

Jo Jn Jo Jn 


m Jo Jn 

because ip vanishes near 912 and near t = T. Finally, 

rT r , r . rT 


[ [ Ue- Pe( [ IVnel^) - [ f f • (f 

Jo Jn ^Jn ' Jo Jn ^Jn ' 


because of (j2.4ip . (I2.45P and the fact that pjz) ^ z for all z > 0 as e \ 0. We thereby see that 
(|1.4p holds and thus infer that u in fact is a weak solution of (ll.ip in 12 x (0, T). The inequality (I2.35P 
results from Lemma 12.61 and the convergence statements. The estimate (I2.36P results from Lemma 
12.71 By (I2.37P and ()2.2p we have the necessary bounds on gradient and initial value, independent of 
E G {Ej)j£^. Furthermore, for any t G [0,T] we can hnd a subsequence {Eji^)k£N of (sjOjeN such that 

and finally the same bound as in Lemma 12.71 holds for u{t) because 

ll«WII$.oo = 


u{t) 

< liminf 

Ue{t) 


oo 

4> 


<liminfmax<^ sup / |Vue(r)|^, ||tto£ - e||<i, oo f 
£=£j>,\0 lo<T<tJn ’ j 


ho||$,oo 


< liminf max < sup / |Vue(r)|^, 

K0<T<tJn 

< max < ess sup / |Vu(r)p, ||uo||<i> c 

I o<T<t Jn 

where for the last inequality we relied on the pointwise a.e. convergence of in (0, T), due to 

()2.45l) valid along a subsequence. □ 

We are now in the position to prove Theorem [la which asserts the existence of a locally positive 
weak solution and Tmax G (0, oo] such that the solution blows up at Tmax or exists globally. 

Proof of Theorem 12.31 By Lemma 12.81 there exists T > 0 such that (jl.ip possesses a locally positive 
weak solution tt on 12 x (0,T), which satisfies (11.71) and (11.81) for a.e. t G (0,T). Hence, the set 


S := <T>0 


there exists a locally positive solution u to (jl.ip on 12 x (0,T) 
satisfying ()1.7I) and ()1.8p for a.e. t G (0, T) 


is not empty and 


Tmax = sup S' G (0, oo] 


is well-dehned. Assume that Tmax < oo and limsup^^j;^^^ ||u(-,t)||^o 


•m 


< oo. 
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This implies the existence of a constant M > 0 such that u < M and hence, due to (HZD, also that 
there is (7 > 0 with < C. Lemma 12.81 provides T > 0 such that for any initial data uq 

satisfying uq < M, < <7, a locally positive weak solution existing on x (0,T) can be 

constructed. 

Choose to G {Tmax — '^,Tmax) such that u{x,to) < M and to)P < C and such that u satishes 

(USD and (jl.Sp at t = to- 

Let V denote the corresponding solution with initial value u{-,to) and define 


u{x,t) = 


u{x, t), a; G il, t < to 

v{x,t — to), X G 0, t G (to, to + T). 


Then u is a solution of m, and (HZD and (11.81) obviously hold for a.e. t G (0, to), whereas for 
t G (to, to + T) we have 

[ |Vu(-,t)|2 

Jn 

< [ |Vu(to)px 

Jn 

( sup [ u{-,t)] ( [ (j)lnu{-,t)- I (j)lnu{-,to) + 
yT-^(^to,to+T)Jn J \Jn' Jn' 

< I iVuopexp ( sup / ^(^r) I ( / (plnu{-,to) - [ (plnuo+ [ 

\Te{o,to)Jn J \Jn' Jn' Jo 


X exp 


rto+T I- 

/ ' 

to J £1^ 

f u 

Jn' 


X exp 


1 


< / iVuol exp 


1 


Jn 


/ (j)lnu{-,t)- / 4>lnu{-,to)+ I I u 

'n' Jn' Jto Jn' 



sup / 


‘JCq' \ T-g(o,t) Jn 


•( / (/>lnri(-,to) - / 4>lnuo+ / u+ / (j)lnu{-,to) + 

^ Jn' Jn' Jo Jn' Jn' Jn' 

= / iVuol^exp 
Jn 



u 


Iq j Q' 



sup / u(-,t) / 4>lnu{-,t)- / 4>lnuo+ I I u 

■^Gfj' \yT&{o,t) Jn J \Jn' Jn' Jo Jn' 

Also, for a.e. t G (0,to + T), 

ll^(•>^)ll<I.,oo < max|||u(-,to)||<j,^^,sup^g(ig_i)^|Vu(-,T)|2| 

< max|max|||no||,j^^,sup^g(o,to)^|Vu(-,r)|2| ,sup^g(iQ^i)^|Vu(-,r)|2| 

< max|||uo||,j,^„o>sup^e(o,b/ol^^(->^)l^} • 

Thus u is dehned on (0,Tmax + ?), contradicting the definition of Tmax- 


□ 


As a direct consequence of (II.8p we obtain that finite-time gradient blow-up cannot occur. More 
precisely, we have the following. 
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Corollary 2.9. Let u and Tmax be as given by Theorem, M.St 
//limsupt/^j.^^^ lk(-,i)llLoo(Q) = oo, then also 

limsup / |Vr(x,= oo. 
t/'Tmax Jn 

Combining now Corollary 12.91 with the estimate (11.71) . we can conclude that if finite-time blow-up 
occurs then also blow-up takes place at the same finite time. 

Corollary 2.10. Let u and Tmax be as given by Theorem \1.3[ 

//limsupt/x^n™, lk(L^)llLo°(o) = oo, then also 

limsup / u{x,t)dx = oo. 
tZ-Tmax Jn 


3. Total mass. Proof of Theorem 11.41 


Let tt be a solution of (HI) on [0, r]. Consider its mass 

y{t) = / u{x,t)dx, t£[0,T), (3.47) 

Jn 

and note that (j3.47jl defines a continuous function on [0,T]. Indeed, we have the following. 

Lemma 3.1. For any weak solution u of a on [0, r], (I3.47P defines an absolutely continuous 
function y. [0,T] —)• M that satisfies 


for almost every t G (0, T). 


y'it) = ivit) - 1 ) / 

Jn 


(3.48) 


Proof. We will show that whenever 0 < s < t < T, 


y{t) - y{s) = {{y{r) - l) JjVu{x, r)|^ dr, 


(3.49) 


where absolute continuity follows from the representation as integral and the assertion about the 
derivative is a direct consequence of division by t — s and passing to the limit s —)• t. 

Let 0 < s < t < T and 0 < d < min {s, T — t}. Define the function y: M —)• M by setting: 

0, T < s — 5, 

1 + s — 6 < T < s, 


X{r) = 1 , 


S < T < t, 


1 - t < T < t + S, 
0, T > t + S. 


Then, according to standard approximation arguments, (p{x,t) := y(t) defines an admissible test 
function for (jl.3p and we obtain 
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Since u G Cioc{[0,T), we have 

rt-\-S 


1 f 1 f 

- / / tt —)• u{t) and - / / n —)• u{s) 

0 Jt Jn ^Js-sJn 


as 6 \0. 

Also by Lebesgue’s dominated convergence theorem, 

pt+S 


and 


rt+5 

Js-S 


[ |Vnp 

Jn 





|Vnp 




as (5 \ 0. Hence, (I3.49p holds. 

This lemma is the main ingredient in the following proof of Theorem 11.41 


□ 


Proof of Theorem \1.4[ (i) In the case of subcritical initial mass Lemma l3. 1 1 shows that y as dehned 
in (I3.47P is decreasing, which by Corollary 12.101 entails global existence, and from the nonnegativity 
of y we derive that y{t) —)• c as t —)• oo for some c > 0. Note that Poincare’s and Holder’s inequalities 
imply that for some Cp > 0 we have 

2 


1 


1 




udx] = 


1 


Cp\n\ 


y 


on (0, oo). 


and hence Lemma [3Tj due to the negativity of y{t) — 1, entails that 


1 


1-2/(0) ..2, 


1 - 2/(0) o 


for almost every t > 0. This would lead to a contradiction to the nonnegativity of y{t) if c were 
positive, whence actually c = 0. 

(ii) If /qUq = 1, then Lemma l3 .1 1 implies that 

rt 


y(i)-l=/ {y{s)-l) [\Vu{x,s)\'^dx 

Jo I Jn 


ds, 


and by virtue of Gronwall’s lemma we conclude y{t) — 1 = 0 throughout the time interval on which 
the solution exists, which combined with Corollarv 1 2.10 1 also implies global existence. 

(iii) In the case when the total mass is supercritical initially, Lemma fS.ll entails that y is nondecreasing, 
and again Poincare’s and Holder’s inequalities imply that 


for a.e. t G [O,^^,^^:) 


^ 2/(0) - 1 2/,N 

with some Cp > 0. Now let z denote the solution to 

z'{t) = ^( 0 ) = ^ 0 , 

for some 1 < 2;o < 2/(0), defined up to its maximal existence time Tq > 0. Then T := Tmax < To-, 
because y > z, and the assertion follows by Theorem 11.31 in combination with Corollarv l2.ini □ 
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4. Global blow-up. Proof of Theorem 11.61 

We proceed to prove that blow-up of our solutions always occurs globally, as stated in Theorem 11.61 

Proof of Theorem \1.6l Assume to the contrary that the closed set ^ is strictly contained in Q. Then 
there exists a smoothly bounded subdomain Q' C such that u is bounded in O' x (0, Tmax)- Let 

4> he a solution to —Acp = 1 in 11', (/> = 0 on cAl'. 

Consider T' < Tmax- Due to the local positivity of u we have ^ € L°°{Tl x (0, T')) and = ^ — 
S x{h,T')) and hence ^ G L2((0, T'), x (0, T')) C L2((0, T'), Wq’^( 11))n 

L°°{Tl X (0, r')). Therefore, it can readily be verified by approximation arguments that it is possible 
to use ^ as a test function in (11.41) . which then leads to 

rt rt r rt 


If- 

Jo Jn' u 


(fdxds + 



Vu ■ V(j)dxds = 


0 Jfi' 


[ ([ (^dx 
Jo \Jn' 


|Vttf dx ds 


for any t G {0,Tmax)- Hence, with Cq' := f^,4> and because of du(t)\g^r < 0, 

/ (l)lnu{t)dx— / (plnuodx— / u ■ Acj)dx ds > Cq/ / / \Vu\'^dxds, 

Jw Jn' Jo Jn' Jo Jn 


that is 


n udxdsp / (^lnu(t)(ix— / (flnuodx > CQ/h{t), 
V Jn’ Jn’ 


(4.1) 


where h{t) := J’^|Vtt(x, s)p dx ds and where - due to the choice of 11' - the left hand side is bounded 
from above. 

On the other hand, from Lemma l3.II we know that 

y'W 


= f \Vu\‘^dx 
Jn 


h{t) = [ / |Vnp dxds = [ f —ds = ln(y(l) — 1) — ln(y(0) — 1) = In 

Jo Jn Jo viV - 1 


t) dx — 1 


y{t) -1 

for y{t) = f^u(x,t)dx. Therefore 
rt r rt 

0 Jn Jo yiv -1 Jn'^odx - 1 

and, by Theorem 11.41 (hi), limsup^^-^^^^ h{t) = oo, contradicting the boundedness of the left hand 
side of (14.ip . □ 

We have seen that the question of global existence versus blow-up of solutions to (II.ip is intimately 
connected with the size of the initial data. If JqUq > 1, the solution blows up globally; if < 1, we 
have proven convergence towards 0. The missing case of solutions emanating from initial data with 
unit mass must exhibit a behaviour different from either, as Theorem O (ii) shows. For a study of 
these solutions, which are actually very important for the described replicator dynamics model, we 
refer the reader to the forthcoming article 


5. Appendix A: Modelling background 

Evolutionary game dynamics is a major part of modern game theory. It was appropriately fostered 
by evolutionary biologists such as W. D. Hamilton and J. Maynard Smith (see [H] for a collection of 
survey papers and [43] for a popularized account) and it actually brought a conceptual revolution to 
the game theory analogous with the one of population dynamics in biology. The resulting population- 
based approach has also found many applications in non-biological fields like economics or learning 
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theory and introduces a significant enrichment of classical game theory which focuses on the concept 
of a rational individual. 

The main subject of evolutionary game dynamics is to explain how a population of players update 
their strategies in the course of a game according to the strategies’ success. This contrasts with 
classical noncooperative game theory that analyzes how rational players will behave through static 
solution concepts such as the Nash Equilibrium (NE) (i.e., a strategy choice for each player whereby 
no individual has a unilateral incentive to change his or her behaviour). 

As Hofbauer and Sigmund m pointed out, strategies with high pay-off will spread within the pop¬ 
ulation through learning, imitation or inheriting processes or even by infection. The pay-offs depend 
on the actions of the co-players, i.e. the frequencies in which the various strategies appear, and since 
these frequencies change according to the pay-offs, a feedback loop appears. The dynamics of this 
feedback loop will determine the long time progress of the game and its investigation is exactly the 
course of evolutionary game theory. 

According to the extensive survey paper [20] there is a variety of different dynamics in evolutionary 
game theory: replicator dynamics, imitation dynamics, best response dynamics, Brown-von Neumann- 
Nash dynamics e.t.c.. However, the dynamics most widely used and studied in the literature on 
evolutionary game theory are the replicator dynamics which was introduced in |50| and baptised in 
m- Such kind of dynamics illustrate the idea that in a dynamic process of evolution a strategy should 
increase in frequency if it is a successful strategy in the sense that individuals playing this strategy 
obtain a higher than average payoff. 

Let us consider a game with m discrete pure strategies, forming the strategy space S = {1,2, 
and corresponding frequencies Pi{t),i = 1,2, ...,m, for any t > 0. (Alternatively S could be consid¬ 
ered as the set of different states (genetic programmes) of a biological population). The frequency 
(probability) vector p{t) = {pi{t),P 2 {t), ...,Pm{t))'^ belongs to the invariant simplex 

S’(m) = jy = (?/i,?/ 2 ,--, 2 /m)^ G K”" : 2 /i > 0,i = 1,2,...,m and ^2/i = l|- 

The game is actually determined by the pay-off matrix A = {aij), which is a real m x m symmetric 
matrix. Pay-off means expected gain, and if an individual plays strategy i against another individual 
following strategy j, then the pay-off to i is defined to be aij while the pay-off to j is aji. For symmetric 
games matrix A is considered to be symmetric. (In the case of a biological population pay-off represents 
fitness, or reproductive success.) 

Then the expected pay-off for an individual playing strategy i can be expressed as 

m 

{A-p{t))i = '^aijPj{t), 
i=i 

whereas the average pay-off over the whole population is given by 

m m 

{p{tf • A ■ p(t)} = aijPi(t}pj(t). 

i=l j=l 

Consider that our game is symmetric with infinitely many players (or that the biological population 
is infinitely big and its generations blend continuously to each other) then we obtain that Pi(t) evolve 
as differentiable functions. Note that the rate of increase of the per capita rate of growth pi/pi of 
strategy (type) i is a measure of its evolutionary success; here pi stands for the time derivative of pi. 
A reasonable assumption, which is also in agreement with the basic tenet of Darwinism, is that the 
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per capita rate of growth (i.e. the logarithmic derivative) Pi/pi is given by the difference between the 
pay-off for strategy (type) i and the average pay-off. This yields the the replicator dynamical system, 


dpi 

dt 


aijPi{t)Pjit) \ Pi{t), i = l,2,...,m, t > 0. 




i=l j=l 


(5.1) 


The dynamical system m actually describes the mechanism that individuals tend to switch to 
strategies that are doing well, or that individuals bear offspring who tend to use the same strategies 
as their parents, and the fitter the individual, the more numerous his offspring. 

Most of the work on replicator dynamics has focused on games that have a finite strategy space, 
thus leading to a dynamical system for the frequencies of the population which is finite dimensional. 
However, interesting applications arise either in biology or economics where the strategy space is not 
finite or, even, not discrete, see [U [35l [36l EZ] . In case the strategy space S is discrete but consisting 
of an infinite number of strategies, e.g. jS = Z, then the replicator dynamics describing the evolution 
of the infinite dimensional vector p{t) = {...,pi{t),p 2 {t),...) is described by the following 


dpi 

dt 


J2aijP,{t)-Y,Yl aijPi{t)pj{t) 

jez i&z 


Pi{t), 


t > 0, 


which is a infinite dynamical system with Pi{t) > 0 for i S Z and ||p(t)|l£i(z) = 1 for any t > 0. 

In the current paper we are concentrating on games whose pure strategies belong to a continuum. For 
instance, this could be the aspiration level of a player or the size of an investment in economics or it 
might arise in situations where the pure strategies correspond to geographical points as in economic 
geography, [26]. On the other hand, in biology such strategies correspond to some continuously varying 
trait such as the sex ratio in a litter or the virulence of an infection, m- There are different ways of 
modelling the evolutionary dynamics in this case, however in the current work we adapt the approach 
introduced in |8|. In that case the strategy set fl is an arbitrary, not necessarily bounded, Borel set 
of M.^,N > 2, hence strategies can be identified by x G 0. For the case of symmetric two-player 
games, the pay-off can be given by a Borel measurable function / : II x —)• M, where /(x, y) is the 
pay-off for player I when she follows strategy x and player 2 plays strategy y. A population is now 
characterized by its state, a probability measure in the measure space (H, £/) where £/ is the Borel 
algebra of subsets of H. The average (mean) pay-off of a sub-population in state against the overall 
population in state ^ is given by the form 

■■= [ [ f{x,y)^{dy)^{dx). 

JQ J fl 

Then, the success (or lack of success) of a strategy x followed by population ^ is provided by the 
difference 

a{x,^):= [ f{x,y)^{dy)- [ [ f{x,y)^{dy)^{dx)=E{5x,^)-E{^,^), 

JQ Jfl Jfl 

where 6x is the unit mass concentrated on the strategy x. 

The evolution in time of the population state ^{t) is given by the replicator dynamics equation 

d£^ r 

——(A)= / a{x, 3{t))^{t)idx), t > 0, cS(0) = 

dt Ja 


(5.2) 
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for any A G j?/, where the time derivative should be understood with respect to the variational norm 
of a subspace of the linear span ^ of The well-posedeness of (I5.2jl as well as relating stability 
issues were investigated in [SGIET] under the assumption that the pay-off function f{x,y) is bounded. 
The abstract form of equation (j5.2jl does not actually allow us to obtain insight on the form of its 
solutions and thus a better understanding of the evolutionary dynamics of the corresponding game. 
In order to have a better overview of the evolutionary game, following the approach in [23t 124] . we 
restrict our attention to measures ^{t) which, for each t > 0, are absolutely continuous with respect 
to the Lebesgue measure, with probability density u{x,t). Then the replicator dynamics equation (j5.2p 
can be reduced to the following integro-differential equation 


du 

'm 


f{x,y)u{y,t)dy- / f{z,y)u{y,t)u{z,t)dydz]u{x,t),t>0, a; G 0, 
1 Jn Jn ) 


(5.3) 


for the density u. 

There are applications both in biology as well as in computer science where the pay-off kernel has the 
form f{x,y) = G{x — y) with G being a steep function of Gaussian type, see [THl [191 [22l [M] . This case, 
in general, models games where the pay-off is measured as the distance from some reference strategy 
and finally under some proper scaling leads to 


/ 

Jn 


fix, y)uiy, t) dy ^ Auix, t), 


(see also [25]) which by virtue of ()5.2p yields 


du 

'm 



u Au dx 


]u. 




(5.4) 


(5.5) 


Another alternative towards getting pay-offs of this type is to consider a game with a discrete strategy 
space and take the appropriate scaling limit. In that case a Taylor expansion and a proper scaling 
gives a similar approximation to (15.4p , see also [23l [24| . 

Therefore in case 14 is a bounded and smooth domain of it is easily seen that via integration 
by parts the nonlocal integro-differential dynamics equation (15.31) is approximated by the degenerate 
nonlocal parabolic equation 

— =u(^Au + J |Vupdx^, X G 14, t>0 

The nonlocal equation (15.6p is associated with initial condition 

u(x,0) = uo(x), X G 14, 

and homogeneous Dirichlet boundary conditions 

u{x,t)=0, X G (911, t > 0, (5-8) 


(5.7) 


(5.6) 


when the agents avoid to play the strategies locating on the boundary of the strategy space since 
they are supposed to be too risky, or the individuals of the biological population do not interact when 
they are close to the spatial boundary where probably the “food” is less. We remark that when on 
the boundary of the strategy space individuals do not really distinguish between nearby strategies 
and hence populate them equally, then the non-local equation (15.6h should rather be complemented 
homogeneous Neumann boundary conditions not explicitly considered here, see |24j . 
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6. Appendix B: A convenient approximation of the initial data 


In the article, we have kept the proof of Lemma l2 .1 1 very short. Here we give a more detailed version, 
which still suppresses some of the more involved technical calculations: 

Proof. Choose 7 > 0 and a domain C ^ such that dist{U^,dQ) > 7 . Let 'd G C^{U^) with i? > 0 
and = 1. Let e > 0 and let cp G be such that \\(p — "WoIIvelzig) < ^ 


$,00 ^ C' + C(e), 

where Q : [0, 00 ) [0, 00 ) is a function satisfying limg^o C(£) = 0- order to see that this is possible, 

recall how smooth approximations p) of VL^’^(n)-functions uq are usually constructed ([58l 13]): With 
the aid of a partition of unity {oi}, the function is written as sum, where the single summands are 
supported in small patches only and those close to the boundary are shifted towards the interior by 
application of shift operators sp, finally the function is smoothened by convolution with a standard 
mollifier . 

We observe that the same procedure applied to <I> does not violate the inequality ||uo||$oo — 
uq < C^, too much, that is 

+ = C$ + C^>, 

holds for some C with hm£\^o C(£) = 0- (The calculations showing this use the fact that mollification of 
smooth functions converge in C^, that $ grows towards the interior, and the Mean Value Theorem.) 
Hence the fact that mollification preserves pointwise estimates that hold everywhere shows that also 
ip satisfies ip{x) < C^{x). 

Let A be a compact subset of Q such that |n \ A| < e and dist(9n, K) < e. Let p G C^{Vl) such that 


p = 1 on A U supp(p and |Vp(x)| < 


£1181(^,90) 


and 0 < p < 1. Denoting 


A = A(e) = [ ^>2|Vp|2 + /(I - p)2|V^>|2 + [ |Vi?|2 ( [{1 - p)$' 

Jn Jn Jn \Jn / 

B = B{e) = - 1 - 2 /(I - p)^> [ V(pVi9 - 2 /(I - p)$ [ {uq - p) [ |Vi9p + 2e|D| /(I - p)$ [ |Vi?|2 

r = r(e) = [ |V<pp + 2 [(uQ-p) [ V<pVi? - 2e|D| [ V<pVi? - 2e\n\ [{uq - p) [ |Vi9p 

J ^ J Q, J ^ J ^ J Q, J ^ 

(Uo-ip)) /|Vl?|2+e2|0|2 /iVlJp 

1. / Jq J q 


+ 


we let C = C{e) = — 


2r 


. Then C solves 


AC^ + HC + r = 0. 


( 6 . 1 ) 


As $ and V<h are bounded, 1 — p is supported on a small set with measure smaller than e, and 
‘hiVpl < 2 D 2 , where $(a;) < D 2 dist{x,dPt), most integrals from the definition of A,B,T can be 
estimated, yielding A —)• 0, H —)• —1, L —)• J(^|VttoP as e ^ 0. Therefore, 

. 2^|Vuo|2 


C = -- 


2r 


= [ |VuoP > 0, 
Jn 


B - \/H2 - 4Ar -1 - 
0, and in particular, limsup(C' — L) < 0. Furthermore, for sufficiently small e, we have C > 0. 


as s 

We also observe that 


a= [{uo-ip)- e|D| -C [{1 - p)^> -)■ 0, 
Jn Jn 
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AS S —0. If e is small enough, therefore, |a| < 


- „„„ .Q - a-RCi 


sup 


hence 


|ai9| < 2 inf{^gf^,dist(x,aD)> 2 } on n (as suppi? C supp(/j). Therefore, 

^p{x) + a'd{x) > - inf ip =: Ck (6-2) 

2 ■(x£f!;dist(x,9f2)>^} 

for X £ K and 

ip + a'd > 0 (6-3) 

on n, because ip > 0 and / 0 only on U^, where (16.21) guarantees (16.211 already. We also have 

p + a'd <{L + ^)<I> + a'd < {L + + d / |tto - 

2 2 

<(L + C(e))‘I> 


with some C fulfilling lim£^oC(£) = 0- Finally, define 

uoe = e + C{1 - p)^ + p{p + ad). (6.4) 

Estimate (j6.3p . the positivity of C and of <I> in fl together with (16.4p entail uqs > £■ Accordingly (12.ip 
holds, for we clearly obtain uoe = e, and Attoe = —C = — /^|VttoeP on dQ, because 

[ iVuoel^ = [ |V(e + C(1 - p)4> + pip + ad)\^ = + {B+ 1)C+ T = C 

Jo, Jo, 

by (|6.1h . Furthermore, 

uoe = / £ + / C{1 - p)^ + pp + a pd= uo 

Jq Jq Jq Jq Jq Jq 

that is ()2.5p . The smoothness assertion follows from the smoothness of p (as mollification) and <I> and 
that oi p,d £ (7^(12). By dehnition of uqs, 

ll^Oe - e||$,oo = l|C'^(l - P)+ Pi^ + «^^)ll#,oo • 

In every point x G 12, uqs — e is a convex combination of and p + ad, which both satisfy the 
estimate “< (L + (C(e))<I>”. Therefore (12.21) holds. Furthermore, 

ll^oe - uo\\wn 2 (^n) = ||e + ^^(l - p) + pip + ad) - uo||tyi, 2 (f^) 

<e^ + C ||V$(1 - p)||^2(f2) + II + C 114.(1 - p)ll^.(^) 

<e\/|f2| + C sup IV$|\/e + 2CD 2 y/s + C sup + e + e + a ||i2||^ 0 

as e \ 0, where we have, once again, used that ||‘hV/3||2^2(Q) < 2 D 2 ^/e, as well as ||'iio|lwi. 2 (f 2 \x) < ^ 
and ||uo — ¥^|IveL2(d) < total, we obtain p2.4p . Finally, given K CC 12, the estimate in (12.3p holds 
for 0 < e < dist(iF, dVl) and with the choice of Ck as in (16.2p . □ 
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